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Dynamical instability of a spin spiral in an interacting Fermi gas as a probe of the 

Stoner transition 
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We propose an experiment to probe ferromagnetic phenomena in an ultracold Fermi gas, while 
alleviating the sensitivity to three-body loss and competing many-body instabilities. The system is 
initialized in a small pitch spin spiral, which becomes unstable in the presence of repulsive interac- 
tions. To linear order the exponentially growing collective modes exhibit critical slowing down close 
to the Stoner transition point. Also, to this order, the dynamics are identical on the paramagnetic 
and ferromagnetic sides of the transition. However, we show that scattering off the exponentially 
growing modes qualitatively alters the collective mode structure. The critical slowing down is 
eliminated and in its place a new unstable branch develops at large wave vectors. Furthermore, 
long-wavelength instabilities are quenched on the paramagnetic side of the transition. We study 
the experimental observation of the instabilities, specifically addressing the trapping geometry and 
how phase-contrast imaging will reveal the emerging domain structure. These probes of the dynam- 
ical phenomena could allow experiments to detect the transition point and distinguish between the 
paramagnetic and ferromagnetic regimes. 

PACS numbers: 03.75.Ss, 71.10.Ca, 67.85.-d 



I. INTRODUCTION 

A magnetic field tuned Feshbach resonance provides 
a powerful tool to control the interaction parameters of 
ultracold atomic Fermi gases [l|. The effective interac- 
tion between two atoms in an s-wave scattering channel 
is attractive on one side of the resonance and repulsive 
on the other side, with both regimes diverging upon ap- 
proaching the resonance. Over the last few years experi- 
ments starting from the attractive side have investigated 
the crossover from a Bardeen-Cooper-Schrieffer state of 
fcrmion pairs to a Bose-Einstein condensate of tightly 
bound molecules [2j. On the other hand, a recent ex- 
periment has provided the first possible evidence for a 
transition to an itinerant ferromagnet beyond a critical 
interaction strength in the repulsive regime [3j. If con- 
firmed, this new realization of ferromagnetism may not 
only resolve long-standing questions stemming from the 
solid state but also promises to open up new arenas of 
ferromagnetism research (j-Q ■ 

When considering the repulsive side of the resonance 
however, it must be noted that the repulsive Fermi gas is 
only a meta-stable state. The two-body ground state in 
this regime is a "Feshbach molecule" , a bound state with 
negative energy. Correspondingly, the many fermion 
ground state is the molecular BEC, whereas the repulsive 
Fermi gas arises only if the system is specially prepared 
without molecules. Even then, atoms gradually recom- 
binc to form Feshbach molecules. At least in the low 
density limit, this occurs predominantly through a three- 
body process @ , whereas closer to the resonance the loss 
may reflect competing many-body instabilities [10| . An 
experiment as in Ref. [3[ must therefore be performed 
inherently out of equilibrium. To abate the fall in atom 
density the experiment was performed while tuning the 



interaction parameter rapidly, but this however masks 
the true phase transition [HI [III- Moreover, it has been 
shown that even if the atom number is kept constant 
by coupling the system to an atom reservoir, the non- 
equilibrium conditions imposed by the three-body loss 
act to change the nature of the ferromagnetic transition 
through an inherently quantum mechanism 13 1. 

In this paper we propose a different strategy to in- 
vestigate the Stoner transition, which could allow inves- 
tigators to circumvent the difficulties imposed by atom 
loss. The idea is to study the dynamical stability of a 
nearly ferromagnetic state, or more precisely a spin spi- 
ral of small wave- vector Q, as shown in Fig. QJb). This 
state has minimal three-body losses as it is locally fully 
polarized. On the other hand, the dynamical stability of 
the spiral spin texture is not protected by spin conser- 
vation because the system has zero net magnetization. 
To take advantage of this new protocol it is first vital 
to determine how the modes of instability change when 
we tune the system across the Stoner transition. For 
an interaction strength tuned so that ferromagnetism is 
favored, the exponentially growing unstable modes are 
expected to reorient the spins, as seen in Fig. QJc), and 
eventually cause the system to fragment into polarized 
domains. A similar collective modes structure was ob- 
served in a bosonic ferromagnetic gas 1141 . Unlike when 
ordering from the paramagnetic state [11], the size of 
these domains and the collective mode structure can be 
finely tuned with the length scale of the initial spin spiral. 

On approaching the Stoner transition we find critical 
slowing down of the unstable modes. However, counter to 
initial heuristic expectations, the unstable modes of the 
helical spin state are, to linear order, the same on the two 
sides of the transition. To differentiate between the fer- 
romagnetic and paramagnetic regimes we go beyond the 
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FIG. 1: (Color online) (a) The gas is started in a fully polar- 
ized state and (b) a normal magnetic field B (purple arrows) 
is applied to form a spin spiral of pitch 1/Q. When the inter- 
action strength is ramped upwards (c) the spins rotate into 
the y direction forming a fully polarized state (which we later 
show has a long wavelength modulation). The thickness of 
the green plane indicates the growing magnetization. In (b) 
and (c) the axis basis set co-rotates with the spin spiral. 



linear analysis and study the feedback effect due to the 
scattering of collective excitations off the exponentially 
growing modes. As the phenomena investigated are not 
only of conceptual interest but could also provide a new 
protocol for the next generation of experiments, we also 
consider the ramifications of a realistic harmonic trap- 
ping potential and the experimental probes of the collec- 
tive modes. 



spin state. The result is a spiral spin texture, as shown 
in Fig. Hlb), with wave vector Q y = (fiB9jt/fi)dB x /dy, 
where gj is the g-factor. The twist rate is independent 
of the spin stiffness and the strength of the repulsive in- 
teractions between particles. This is in close analogy to 
the effect of a potential gradient placed across a super- 
fluid, which is essentially a XY ("phase") ferromagnct. 
The potential gradient affects phase twist at a constant 
rate, independent of the supcrffuid stiffness, which cor- 
responds to free acceleration according to Newton's law. 

To gain further understanding into why interactions 
do not impact on the dynamic formation of the initial 
spin spiral we can study the system within the simple 
setting of the Heisenberg ferromagnet. Then, the situ- 
ation posed in Fig. Q] is described by the Hamiltonian 
H = -Jj2{ij) Sj • Sj - bYsiViSf, where the first sum- 
mation covers nearest neighbors on the lattice, J denotes 
the coupling between adjacent sites, and b is the magnetic 
field gradient. The initial conditions are Sf(t = 0) = 0, 
SJ (t = 0) = 0, and Sf(t = 0) = S . We then study 
the evolution of the spins using ftS = i[if, S], finding 
the equations of motion hS x ^ = 0, HSy t i = byiS z j, and 
fiS x ,i = —bytSy.i, where the components containing J 
cancel exactly. Finally, we can then solve the equations 
of motion to yield S Xt i(t) = 0, S y .i(t) = Sosin.(byit/H), 
and S Zt i(t) = So cosibyit / K) . These show that the for- 
mation of the spin spiral by the external magnetic field 
gradient is independent of the interactions, J, between 
particles. 

An interesting and beneficial practical implication of 
the above observation is to alleviate the need to first form 
a spiral, and then perform a Feshbach field quench. The 
spin spiral may equally well be formed with the Feshbach 
field in place. Finally we note that a magnetic field gra- 
dient has the side effect of imparting a translational force 
on the gas perpendicular to the gradient. Fortunately, in 
this scheme we need to impose only a long-pitched spin 
spiral, which requires only a weak magnetic field gradient 
with minimal side effect. 

Following the preparation stage, we arc ready to allow 
the spiral state to evolve under the influence of the repul- 
sive interactions tuned by the Feshbach field, and track 
its evolution into the polarized state out of the plane of 
the initial spiral shown in Fig. [He). 



II. FORMING THE SPIN SPIRAL 

To form the initial spin spiral, the atomic gas is first 
prepared in a fully polarized phase, say along the z spin 
axis shown in Fig. [TJa), and a magnetic field gradient 
is imposed perpendicular to the magnetization axis (e.g. 
B = byx) for time t. Such a field can be thought of 
as a gradient of the relative potential between the fx 
and 4-x spins, and leads to a (constant) relative accelera- 
tion between these two spin components in the coherent 



III. LINEAR SPIN- WAVE INSTABILITY 

We now investigate the instabilities of the helical spin 
state by first focusing on the linearized spin fluctuations 
around the initial spin spiral. The exponentially growing 
unstable modes will show up in this analysis as collective 
excitations with imaginary frequencies. As the dominant 
modes grow exponentially, non-linear processes become 
important. Later in Sec. IIVI we will study how scatter- 
ing off the exponentially growing modes renormalizes the 
spectrum. 



To study the collective modes we start from the quan- 
tum partition function expressed as a fermionic coherent 



state path integral, Z = Tr C -P( H -^ N ) 
the corresponding action 
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where J = f„ dr J dr, the free particle dispersion is 
£k = fc 2 /2, /J. is the chemical potential, and <?<5 3 (r) is 
the strength of the s-wave repulsive contact interaction. 
We have also set H = m = 1. To explore how interac- 
tions impinge on the collective mode spectrum we affect 
a Hubbard-Stratonovich decoupling, which incorporates 
the spin channels cf> + <p 0, [l5( ■ The static component 
of the magnetization, </>o, follows the initial helical spa- 
tial spin texture and <p represents the growing unstable 
modes and fluctuations around that stationary compo- 
nent. It is also convenient to apply a gauge transforma- 
tion ip i— > ^c-'Q' 1 ""*/ 2 to enter a spatially rotating basis 
set with pitch vector Qe x /2, which renders the initial 
spin texture, and thus also the static component of the 
magnetization to be uniform, 0o = (0, 0, 4>o). 

After integrating out the Grassmann fields we obtain 
Z = j V<pc^ s with the action 



S=|.g^-Trln(G- 1 ) 



g<p — Trln 
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where Gq = d T + e^ — fi — ga z (f)o = G~^ denotes the 
elements of the inverse Green function at the level of the 
renormalizcd mean field. We then expand in the pertur- 
bative collective modes <fi. We assume that the initial 
spin spiral has a long wavelength relative to the Fermi 
wave vector fcp, so we work in the regime Q <C kp. Fo- 
cusing on the soft modes that are perpendicular to the 
saddle point field, the dispersion satisfies u) -C \i and has 
wave vector q -C kp. Following the expansion for the 
soft modes we get the contribution to the action from 
fluctuations 
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Before searching for the collective modes it is useful 
to transform the basis set for the magnetization from 
Cartesian to spherical polar coordinates: <f> = 0o[l + 
rj\(cos 6, sin 9 sin tp, sin 9 cos tp). Note that 6 is defined as 
the angle with respect to the positive x axis, and tp is 
the angle from the z axis of the projection of 4> onto the 



y,z plane. The initial state is 9 = 7r/2 and tp = 0, which 
ensures that fluctuations in <p are well defined. With this 
definition we expand in small thermal and quantum fluc- 
tuations away from the fully polarized state that are both 
angular tp and 6 (now redefined to be centered around 
7r/2), and also in the amplitude of the mode, r\. The per- 
turbative form for the expansion correct to quadratic or- 
der is <p = 4> Q [l+rj\{-e,tp, l-[6» 2 + v? 2 ]/2). Although both 
Cartesian and spherical polar basis sets yield the same 
collective mode structure in this linear response analysis, 
when in Sec. II VI we consider feedback corrections to this 
response it will be necessary to work with the spherical 
polar basis set to properly evaluate phase and amplitude 
fluctuations. However, to study the system with only lin- 
ear response we also expand G± in uj and q and find that 
the soft modes are coupled in the action through 



S = Sn Tr 
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!df> i 2 = i 2 



) + g 2 , and 5 = 

2 1 / 3 /3(/>ofcFa/7r. The twist wave vector Q couples only 
to the magnetization along the initial spin spiral. Due 
to their commutation relations, the conjugate modes are 
coupled by the off-diagonal elements and so to extract 
the dispersion we demand that the determinant of the 
matrix is zero. This yields the dispersion 



±xqV<i 2 - Q 2 



(5) 



which is shown in Fig. [21 a). We first verify the dispersion 
in the absence of the spin spiral. When Q = and in- 
teractions are strong so %-> 1/2 we recover the familiar 
dispersion of a single particle, ui = q 2 /2. When q < Q 
the dispersion is imaginary, corresponding to an instabil- 
ity, whereas when q > Q we recover oscillating modes. 
Note that, as required by spin conservation, there is no 
dynamical instability of the magnetization at zero wave- 
vector. Similarly growth at q — Q is stunted as this 
mode is initially fully polarized. The exponential growth 
of the order parameter is maximal at the wave-vector 
q = Q/y2, and it is at this length scale that we would 
expect to see ferromagnetic domains emerge. In the ex- 
periment [3[ domain walls could not be observed in a 
gas starting from a paramagnetic state since their size 
falls below the resolution of current experimental tech- 
niques 11 1, but here their length scale can be tuned to 
be experimentally observable by changing the pitch of the 
spin spiral. The fully polarized phase becomes unstable 
at the critical interaction strength kpa = 2 5 / 3 3/5 w 1.90. 
At this interaction strength x = so the system under- 
goes critical slowing and the characteristic time of the 
instability diverges. 

Against intuitive expectations we find the same collec- 
tive mode structure, to linear order, on either side of the 
transition. To gain insight on this, it is useful to study 
the linear stability of a spiral texture in the Heisenberg 
model on a lattice H = —JY^Uj) Si ■ Sj. Here (ij) re- 
stricts the summation to nearest neighbors and J denotes 



the coupling between adjacent sites. A straightforward 
spin wave analysis of the fluctuations about the spiral 
texture yields the collective mode dispersion 



cosi(Qa) sin 2 (^)+sin 2 (^)+sin 2 (^ 



x f m {—) +sm {—) +sm {—)- sm {—) ' (6) 

where a is the lattice spacing. Note that in the limit 
qa <C 1 and Qa <C 1 the collective mode dispersion 
Eq. ^) computed in the Heisenberg model approaches 
the same form Eq. ([5]) developed in the continuum case. 
Tuning the coupling from ferromagnetic to antiferromag- 
nctic through criticality at J = is allied with a vanish- 
ing dis[persion, as found at the Stoner transition of the 
Fermi system. The anisotropy introduced by the first 
root gives a preference to modes with a wave vector in 
the plane of the spiral. 

We note that the equivalent dynamics of the initial 
spin-spiral seen for positive and negative J is due to an 
exact symmetry of the Heisenberg Hamiltonian. This can 
be shown by a variation on the arguments presented in 
Ref. [16[, which is valid for initial states invariant under 
some generalized time reversal transformation. Such an 
exact symmetry does not exist in the itinerant fermion 
model of interest here. Rather the symmetry of the col- 
lective mode dispersion around the dynamical critical 
point is an emergent phenomenon, and as we shall see 
later in Sec. IIV1 valid only within the linear approxima- 
tion of the dynamics. 

Before proceeding to study these nonlinear effects, 
we revisit our initial assumption that the three-body 
loss may be neglected in the proposed setup. To as- 
sess the validity of this assumption we shall compare 
the loss rate in the spin spiral to the growth rate of 
the maximally unstable mode around this spiral state. 
The three-body loss rate (strictly valid for fc F a <C 1) 
is llle(fc F a) 6 n-j-n^(n-j- + n\) [||. To evaluate this we 
note that adjacent spins in the spiral, separated by 
r s = (2/97r) 1 / 3 fcp 1 , have an angle between them of 
Qr B -C 1. Therefore the adjacent spins give a geomet- 
ric component of n^n± = (2/97r) 2 / 3 n 2 (Q/fc F ) 2 /4. Apply- 
ing this to the experimental regime gives the loss rate 
0.12(fc F a) 6 (g/fc F ) 2 s- 1 . 

On the other hand, the growth rate of the domi- 
nant mode, according to Eq. ([5]), is max q {Im[w(q)]} = 
(1 - 2 5 / 3 3/5fc F a)Q 2 /4. Interestingly we observe that the 
ratio of loss to growth rate depends only on the inter- 
action strength and not on the spiral pitch Q. Com- 
paring the two rates we find that the loss will dominate 
for strong interactions fc F a > 4.9, though in this regime 
the formula for loss rate significantly overestimates the 
true loss [13] and in fact the theory should be valid for 
even higher interaction strengths. However, in the exper- 
imentally accessible region including the phase transition, 



k-pa < 2.5, loss is more than 60 times smaller than the 
dominant growth rate of the collective modes. Therefore 
the new experimental protocol offers a promising way to 
observe magnetism without the damaging effects of loss. 
In addition, losses will dominate within a small region 
of width d(k-pa) ~ 0.03 surrounding the point of critical 
slowing down. However, we shall sec in the next section 
that non-linear feedback effects drive a new dynamical 
instability at large wave- vectors precisely in the region of 
critical slowing down, and that this instability will over- 
come the losses. 



IV. NONLINEAR COLLECTIVE MODES 

In the previous section we found that the behavior of 
the collective modes spectrum with interaction strength 
is qualitatively similar cither side of the critical interac- 
tion strength fc F a = 2 5 / 3 3/5. However, for the Fermi sys- 
tem this is only a feature of the linear analysis, which we 
now extend to include, self-consistently, the effect of the 
scattering of fluctuations on the exponentially growing 
modes. We will see that this mechanism also eliminates 
the region of critical slowing down by generating new 
unstable modes at high wave- vectors near to the critical 
interaction strength. 

We focus solely on the consequences of scattering off 
a single dominant growth mode 3?p. According to the 
analysis of Sec. Hill the momentum P associated with the 
dominant mode is Qj\f2. though here we allow a general 
wave vector that we will later determine self consistently. 
In presence of the large mode $p , we have to expand the 
action Eq. ([2]) to cubic order to include terms which are 
linear in $p and quadratic in the other modes 

S* = f g 3 Tr [G_ ($ x + i$ y )G+</> z G+(0 x - i<f> y ) 
- G + ($ x -r$ y )G_0 z G_(</> x + i0 y )] , (7) 

where G± = 9 r + eg.— /i=pg0o- This new scattering mech- 
anism couples the {(/i x , </> y } channels to the <j) z channel, 
and so describes scattering out of the original magnetiza- 
tion configuration and into the dominant growing mode. 
The presence of a growing classical field $(£) = $(0)e iSlr 
(with imaginary frequency 0), requires a seed fluctuation 
in the initial spiral state. Such a seed will be present in 
any realistic implementation because of random inhomo- 
geneity in the magnetic field and thermal fluctuations. 

We again transform to the spherical polar basis set 
to properly separate phase and amplitude fluctuations. 
Coupling to the exponentially growing mode of wave vec- 
tor P and frequency SI requires studying the spin suscep- 
tibility matrix expanded out to include couplings between 
other modes at wave vectors q and q±P, and frequencies 
w and ui ±Q. The coupled action then takes the form 
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FIG. 2: (Color online) The variation of the exponential growth rate of the collective mode at wave vector q with interaction 
strength kpa. Calculated for an initial seed field of size (a) $ = 0, (b) $ = 0.1Q 2 /g, and (c) $ = 0.2Q 2 /g. The fields {9, ip} 
are evaluated at wave vector q, and r\ at wave vector q± Q; we focus on the peak contribution to the growth with q || Q. The 
black trajectory highlights the maximum growth rate at a given interaction strength. 
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where £ = 3Ep8 s (l — ^ fc ^ )/2, -Ep is the polarized state 
Fermi energy, and 9 S denotes the seed field size for ei- 
ther 9 or tp. The expansion in frequencies employed for 
the amplitude modes, 77, applies for w <C q, which holds 
for the regime of interest that describes the coupling the 
soft angular modes to the amplitude modes. We again 
search for the zeros in the resultant determinant to ex- 
tract the collective modes and explore the unstable re- 
gion. Note that in the absence of the growing field, the 
formalism immediately recovers the dispersion found in 
Sec. IIHI However, here we aim to go further and consider 
the form of the modes in the presence of the dominant 
growing field. To ensure that the formalism is self con- 
sistent at each interaction strength we first compute the 
feedback-corrected dominant mode before calculating the 
entire collective modes spectrum in the presence of that 
growing mode. We also note that this formalism respects 
the the spin conservation law, and so, as in the linear 
analysis of Sec. IIII1 growth of the uniform component of 
the magnetization is suppressed. 

The final result for the spectrum is shown in Fig. [2] for 
different values of the seed field. First we note that for 
interaction strengths far from the critical slowing, scat- 
tering on the growing mode is negligible. In particular 
the dominant growth mode remains at ~ Q / V2 in this 
regime. However, near to critical slowing the wave vector 
of the dominant growth mode can be enhanced. This can 



be understood by comparing the on and off-diagonal ele- 
ments of the action matrix, Eq. (|8J) . Upon nearing criti- 
cal slowing, the on-diagonal elements, being proportional 
to Xj approach zero. However, the off-diagonal terms, £, 
that represent scattering off the dominant growing modes 
are non-zero. To ensure that the overall determinant is 
zero demands a non-zero value for the on-diagonal ele- 
ments which requires a large wave vector. On approach- 
ing critical slowing this growth in the wave vector would 
only be curtailed by higher order momentum terms. A 
detailed analysis shows that the wave vector of the peak 
growth tracks the trajectory along the critical slowing at 
kpa = 2 5 ' 3 3/5 in the phase diagram Fig. [2] The peak 
growth at a particular wave vector q is when that wave 
vector q is parallel to the direction of Q, and so we focus 
on that contribution in Fig. [2] On the paramagnetic side 
of the resonance, kpa < 2 5 / 3 3/5, the enhanced scatter- 
ing into the higher momentum sector removes long wave- 
length components of the domains, whereas on the ferro- 
magnetic side of the resonance, k-pa > 2 5 / 3 3/5, the larger 
domains are naturally still favored. Therefore, consider- 
ation of the feedback of the exponentially growing field 
yields an additional collective mode structure that could 
experimentally distinguish between the two sides of the 
critical slowing interaction strength. 

This picture was developed in the presence of just a sin- 
gle growing exponential mode, whereas in reality these 
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FIG. 3: (Color online) The variation of growth rate (red solid, 
primary y-axis) and domain size (blue dashed, secondary y- 
axis) with radius within a harmonic well, (a) corresponds 
to linear response, and when (b) feedback corrections with a 
seed field of size 0.2Q 2 /g are taken into account. The vertical 
cyan line denotes the radius corresponding to critical slowing 
and r m ax is the outer radius of the atom distribution. 



seed modes have an initial growth rate ~ XQvQ 2 ~ Q 2 
(Eq. ((5)) centered around the wave vector q = Q/V2 
corresponding to maximal growth. With reference to 
Fig. [2J except near to critical slowing, at a given inter- 
action strength the mode growth rate is sharply peaked 
as a function of q so that the presence of other less dom- 
inant exponentially growing modes blurs the collective 
modes dispersion by less than 5%. The relative initial 
size of those modes can however have an impact on the 
mode spectrum. Created by inhomogeneities in the ap- 
plied magnetic field, the uncertainty can correspond to 
seeing a range of growth rates, for example the range 
spanning between Fig. HJa) and Fig. |5fc). Though the 
growth rate around critical slowing is quantitatively dif- 
ferent, the qualitative features are robust and the dom- 
inant growth wave vectors remain the same, so in the 
experiment strong signatures of ferromagnetism should 
be observed. 



V. EXPERIMENTAL OBSERVATION 

Having studied the instability to a ferromagnetic state 
in a uniform gas we now turn to consider the experimen- 
tal ramifications of our results. Wc will focus on two key 
questions: first the consequences of the atomic gas be- 
ing held within a realistic harmonic trapping potential, 
and second we address the experimental signatures of our 
predictions. 

Wc shall treat the harmonic confinement within the lo- 
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FIG. 4: (Color online) The Fourier transform of the domain 
distribution within a harmonic well. Without magnetization 
feedback is shown by the red solid line, and the blue dashed 
curve is in the presence of magnetization feedback. The peak 
growth rate at Q/y/2 is labelled by the vertical cyan line, 
whereas the maximum growth wave vector Q is shown by the 
vertical green line. 



cal density approximation. On passing radially outwards 
from the center of the harmonic well the local density and 
therefore the effective interaction strength falls. We can 
therefore use the results of the previous sections to map 
the rate at which the instability develops and its char- 
acteristic wavelength as a function of the radius. These 
results are shown in Fig. [3] for both the linear analy- 
sis (a) and when the feedback due to scattering on the 
growing mode is taken into account (b) . Note significant 
modifications due to feedback corrections. First near the 
radius which corresponds to critical interactions there is 
a distinct maximum rather than a vanishing growth rate. 
The enhanced growth rate is driven by a scattering off 
the growing mode into large wave vectors. Therefore, 
secondly the characteristic wavelength of the unstable 
modes dips around that radius. In the linear analysis, 
by contrast the characteristic wavelength of the unstable 
modes is everywhere 2n\/2/Q. 

We now turn to the question of how the distinct spatial 
structure of the instability may be observed in an exper- 
iment. The first approach we look at is to employ dif- 
ferential in-situ phase-contrast imaging, a method that 
has already been used on the ferromagnetic cold atom 
gas Q . This measures the difference between the up and 
down-spin populations integrated along vertical columns 
through the gas. If the procedure is repeated across the 
gas, a two-dimensional map will be produced, which is 
governed by the pattern of the emergent magnetization 
aligned normal to the imaging plane. Once Fourier trans- 
formed, the spectrum of this map will reveal the typical 
magnetic domain size. To simulate the expected exper- 
imental result we employ a heuristic model of the gas. 
We divide the system into domains, with position de- 
pendent size, determined by the wavelength of the dom- 
inant growth mode at that radius. Each domain is then 
assigned cither up or down magnetization at random, 



and the resulting magnetization structure is column inte- 
grated to obtain a two dimensional magnetization image 
as in the experiment. The power spectrum of the do- 
main structure in the simulated image is shown in Fig. 3) 
Without feedback corrections, according to Eq. ([5]) the 
dominant mode is at q = Q/\/2 irrespective of inter- 
action strength. In the numerical experiment a distinct 
peak exists at wave vectors around and below Q/y2; this 
is because adjacent domains were not necessarily mis- 
aligned thus increasing the effective length scale of the 
observed domains. If feedback corrections are taken into 
account then in the region of the trap corresponding to 
critical slowing significantly smaller domains are formed. 
This is reflected in the Fourier spectrum with reduced 
weight at small wave vectors, and enhanced weight at 
larger wave vectors. By tuning with cither the size of the 
seed field or the duration of the experiment, the signal of 
the contrasting behavior when feedback corrections are 
taken into account could help investigators to identify 
the ferromagnetic transition. 

One criticism of the recent experiment [3| which re- 
ported the first signs of ferromagnetism is that the do- 
mains, if present, were too small to image. In the new 
experimental protocol the pitch length of the initial spin 
spiral sets the size of the domains formed. However, if 
this length scale remains below the resolution of the ex- 
periment then a statistical analysis of the in-situ phase- 
contrast imaging could still provide an estimate of the do- 
main size. If the orientation of adjacent domains within 
a column is independent, then measurements over ad- 
jacent columns will give an estimate of the variance of 
magnetization in each column. The uncertainty in the 
net magnetization should vary as sN/A^/n, where s is 
the spin per atom, N the number of atoms imaged, and 
n the number of domains. A larger number of domains 
will result in a smaller variance in the net magnetiza- 
tion, thus allowing the total number of domains to be 
estimated. 

One final experimental probe is spin-dependent Bragg 
spectroscopy. This could track the decay of the planar 
spin spiral and the reducing signature at wave vector P. 
Furthermore, using a variable wavelength optical lattice 
potential to couple asymmetrically to the spin degrees of 
freedom, the collective mode response could be studied 
through dynamical fluctuations of the cloud spatial dis- 



tribution as a function of wavelength, laser amplitude, 
and detuning. 



VI. CONCLUSIONS 

We have studied the time evolution of a spiraling spin 
texture, prepared by a magnetic field gradient in an in- 
teracting degenerate Fermi gas. The linearized dynamics 
of the magnetization shows an instability, which devel- 
ops canting out of the spiral plane, with the maximally 
unstable mode at a wavevector l/v2 that of the origi- 
nal spiral. The instability grows exponentially in time 
with a characteristic time scale that vanishes close to the 
Stoncr transition point. To linear order, the dynamics is 
the same on the ferromagnetic and paramagnetic sides of 
the critical slowing down. 

Interestingly, however, we find that near the critical 
point, non linear effects in the dynamics become impor- 
tant. Specifically, scattering of collective excitations on 
the exponentially growing unstable mode acts to renor- 
malize the spectrum, shifting the instability to larger 
wave- vectors on approaching the original critical point. 
Moreover the point of critical slowing is eliminated and 
the new branch of excitations at high wave-vectors al- 
lows a clear distinction between the ferromagnetic and 
paramegnetic regimes. 

Finally, we studied the experimental signatures of the 
unstable dynamical modes in a realistic trap confinement. 
Using the calculated dispersion of unstable modes we ob- 
tained the spatial distribution of spin-domain sizes in the 
trap, which may be observed with phase contrast imag- 
ing. The most dramatic signature of the new branch of 
instabilities induced by the non linear feedback effect is a 
collapse of the domain size at a particular critical radius 
in the trap. 
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